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1 Introduction 

O 

In many areas of practical interest, e.g. computational fluid dynamics or image reconstruc- 
tion, the computations are complicated and expensive, effectively limiting the achievable 
precision. In order to overcome these problems, adaptive finite element approaches are in 
use since several decades (see, e.g., [BR78 ). For instance, in the so-called /i-adaptive meth- 
ods the computational meshes are refined locally so that the mesh captures the variation 
of the solution while remaining coarse elsewhere. It has been shown that such approaches 
are computationally much more efficient than uniform meshes. In recent years, there has 
^ \ been considerable progress in applying these techniques to more involved questions such as 

the a posteriori error estimation of values of nonlinear functionals of interest {goal- oriented 
estimation, see, e.g., |BR01| . [BR03| . [Ran05| ) or the (additional) a posteriori estimation of 
modeling errors (see, e.g., |OV00j . [VOOlj . |BE03] ). 

The present paper describes the extension of recent techniques for obtaining a posteriori 
error estimates for modeling and discretization errors to nonlinear second-order elliptic PDEs 
which are discretized by means of node-centered finite volume schemes including stabilization 
mechanisms of upwind type. Finite volume methods are attractive methods in selected areas 
of application, and therefore it is a natural requirement to develop analogous methods of error 
control for FVM. However, since finite volume methods suffer, in general, from the so-called 
property of Galerkin-orthogonality, special attention is to be paid to the treatment of the 
resulting defect term. It is shown that the extension of the dual-weighted a posteriori error 
estimates to finite volume discretizations is possible in a reasonable way. Furthermore, the 
latter approach is interesting because of the fact that different sources of errors (i.e. not only 
discretization errors but, for example, also modeling errors) can be estimated with respect to a 
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rather arbitrary user-defined output functional. For instance, in the field of inverse problems, 
the Tikhonov functionals can serve as typical output functionals (see, e.g., jBKlO]). 

Here we will mainly deal with Voronoi and Donald finite volume partitions on simpli- 
cial primary partitions of the computational domain; however the ideas can be extended to 
more general primary partitions, in particular quadrilateral or hexahedral partitions (cf., e.g., 
|Ang06j Sect. 4.2]). 



We consider the following boundary value problem with respect to the unknown function 
u: Q^R: 

-V- (A(-,u)Vu) + b(-,u) • Vu + c(-,u)u = f in Q, , . 

a = onT, U 

where f2 C R d , d £ {2, 3}, is a bounded polygonal or polyhedral domain with a Lipschitzian 
boundary T, and the data in ([I]) are sufficiently smooth: 

A : fixK-)> R d ' d , h: nxR->M. d , c: (Ixl^M, / : fi R. 

Equations of such type may occur in various areas of science, for example in the mathematical 
description of filtration processes in nonhomogeneous media. 
Using the formal notation 



(w,v) := / wvdx, 
Jn 

(Vu>,Vt>) := / Vw-Vvdx, 
Jn 



b(w;v) := -[(b{-,w)-Vw,v)-{wM-,w)-Vv)), (2) 
d(w; v) := (c(-,w)w, v) - - (u>V • h(-,w),v) , 

a(w;v) := (A(-, w)Vw, Vv) + b(w;v) + d(w;v), (3) 

and (/, v) := {f,v), the variational formulation of the problem ([T]) in the space V := Hq(U) 
reads as follows: 

Find u G V such that 

ViiGl/: a{u;v) = (f,v). (4) 

Regarding results for the existence, uniqueness and regularity of so lutions of ([1]) or (JU), there 



is a wide literature both of relatively general nature (see, e.g., [B08[ Ch. 2] for a short survey) 
as well as for more specialized equations (see, e.g., |AS06| ). 



2 The finite volume scheme 

Finite volume methods are attractive discretization methods for partial differential equations 
of first or second order in conservative form since they adequately transfer the conservation 
law, which is expressed by the differential equation, to the discrete level. At the same time, 
due to their proximity to finite difference methods, they are relatively easy to implement even 
in the nonlinear situation. However, a drawback of many finite volume methods is that there 
is no p- hierarchy as in finite element methods, therefore the order of accuracy (related to the 
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Figure 1: Configuration for the Voronoi-type discretization (d = 2) 



grid size) is relatively low. Nevertheless finite volume methods find wide applications in the 
computational practice. A certain degree of compensation can be achieved by the application 
of adaptive techniques based on a posteriori error estimates, as discussed in the subsequent 
section. 

In this section we concentrate on node-centered finite volume methods for the discretiza- 
tion of problem (pTJ) . 

2.1 The case of Voronoi diagrams and scalar diffusion coefficients 

Let us consider a family of Voronoi diagrams such that their straight-line duals are Delaunay 
triangulations of Q consisting of self-centered simplices. Here a simplex T is called self-centered 
if its circumcentre lies in the interior of T or on the boundary dT. 

Denote by A C N the index set of all vertices Xi of a particular triangulation T and by 
A C A the index set of all vertices lying in Q. 

In more detail, let 



:= n* 




{ie(l: \\x — Xi\\ < \\x — Xi\\ V j G A \ {i}}, i G A, 
where || • || denotes the Euclidean norm in M. d , 
measd (fi*) , 

where meas^ (•) denotes the d-dimensional volume, 

dtti n dn j} rjj := n t, % e A,j e A \ {i}, t g r, 

meas rf _i (Iy) , mfj := meas rf _i (r^) , 

1 1 1 1 j 



h 



{j eA\{i} : m^O}, 

{i G A : G ar}, 

max/17 1 , where /it := diamT. 
TeT 
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The finite volume solution will be interpolated in the discrete space 

V r :={v€V: (VT G T : u | T G Pi( T ))} , 

where V\ (T) is the set of all first degree polynomials on T. We introduce a so called lumping 
operator 

Lj : C(Q) ->■ Loo(^) acting as := '^2v(x i )xn i , 

ieA 

where Xfii denotes the indicator function of the set fij. 

Due to stability reasons, especially for the case of dominating convection, the class of finite 
volume methods under consideration is characterized by an additional stabilization technique 
called upwinding. For that purpose we introduce a scaling function K : R — > [0, oo) which is 
defined by the help of a weighting function r : R — > [0, 1] as K{z) := 1 — [1 — r(z)]z. 

A typical example of such a weighting function is 

leading to K(z) = z/{e z — 1), the Bernoulli function. 

The discrete problem for the case of a scalar diffusion coefficient, i.e. where A is of the 
form AI with A: x R — > R and I being the identity in R d , is formulated as follows: 

Find uj- € Vqr such that 

VV e : ar(«r; «r) = (/r, «r), (6) 




(u>Ti - u>Tj)mij + CiWrifrii 



Xj + Xj 107-j + Wfj 

2 ' 2 



lij = Jij{w T i, w T j) := Vij ■ b I , 1 

Ci = Ci(w T i) := c(xi,w T i), fi-=f(xi). 



Moreover, we introduce the following norms and seminorms on Vj- '■ 

\\vt\\t ■= V {vr, vt)t = 1 1 L T vr 1 1 o,2,n, 



1/2 



i n r i 1 2 i ii ii2 1 V 2 

\vr\\v ■= {WW + \\ v r\\r\ ■ 
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For the sake of consistency in the notations, we also use the following abbreviations of well- 
known seminorms and norms in the Sobolev space H 1 ^): 



\vt\d '■= \VT\l,2,n, 



\ v t\\d 



\VT\\i,2,n- 



The scheme © with the weighting function ([S} is often called exponentially upwinded. It can 
be defined for other control functions r : R — > [0,1], too. However, we have to assume that 
all of these control functions satisfy the following properties: 



(PI) 
(P2) 

(P3) 
(P4) 

(P5) 



r(z) is monotone for all z G R, 
lim r(z) = 0, lim r(z) = 1, 

z— >■— oo z— >oo 

1 + zr(z) > for all z G R, 

[1 — r(z) — r(— z)]z = for all z G 



r(z) 



z > for all z G 



(P6) zr(z) is Lipschitz-continuous for all z G R. 
We get from (1P4|) the relation 



l + zr(z) =K(-z). 
Replacing in ([7|) the argument z by — (JP3J) immediately implies 
(P7) K(z) > for all z G R. 

EXAMPLE 1 T/ie function 



(7) 



-[signz + 1], 



due to JBT81J, has been investigated in \Ris86\j , lRis90f for a linear equation in £7]). This 
scheme is called fully upwinded. 

The next two examples are simple approximations of ©. 
EXAMPLE 2 



EXAMPLE 3 



r(z) 



r(z) 




-m 
m 
rn 



< m < 8, 



r o , 


z 


< 


—m 


0.5 , 


\z\ 


< 


m 


> 1 . 


z 


> 


m 


tm. 









< m < 2. 



EXAMPLE 4 The choice of the function 

r(z) = - 

goes back to Samarskij [Sam65]. 



2 + \z\ 



+ 1 
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EXAMPLE 5 In \MCL+ W$ the function 

r(z) = — [tanh z + 1] 

was proposed. 

EXAMPLE 6 Finally, it can be taken the function 

[1-<t)/2 , z<0 



r(z) 



with a(z) := max <j 0, 1 — y- - f , 



[l + a)/2 , z>0 
what corresponds to Ikeda's partial upwind scheme E JTke83]. 

Both for computational and theoretical reasons it may be advisable, in a really nonlinear 
situation, to choose differentiable control functions r. In the sequel, if there is no special 
reference, we assume that the scheme under consideration is defined for a general function r 
that possesses the properties ()Pip to (|P6p . 

Finally we mention two equivalent representations of the form a-j. First we remember 

that the leading coefficient ^p~K ' 

dij 



the following manner: 



where r^j := r 



flij 
lijdij 



flij 
lijdjj 



in aj- can be written, by the definition of K, in 

lij dij 



1 — r 



f'ij 



d, , 



(l-rij)^ 



. Hence we get the representation 



fHj 

a T {w T ; vr) 



ieA 



dij 



(1 - Tij) (w T i - "'7 , )"-,'"/./ } 



+ CiW Ti m^. 
Furthermore, introducing the notations 
a? r {w T \ v T ) ■ 



(8) 



ieA jeAi 



b r (w r ; vr) 
d T {w T ; vr) 



■ieA jeAi 




ami 



(l - nj)w T j 



\ E ^i mi i \ 
jeAi ) 



Tij W T i 



Tij TTT'ij j 



WTiVTi 



we get a splitting of ar which is comparable with (|3j): 

a T {w T ; vr) = a T (w T ; vr) + M™r; vr) + d T (w T ; v T ). 



(9) 
(10) 

(11) 
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REMARK 1 In the special case V • b = on Q, it is senseful to use the following versions 
of b-r and d-r: 



b T (w T ; v T ) = y^VTi ^ [(1 ~ r ij) w Tj + rijW Ti ] jijirii 
ieA jeAi 

d T (w T ;vr) = } ] CjWriVTi'mi- 
ieA 

2.2 The case of Voronoi diagrams and matrix-valued diffusion coefficients 

Using the representation 

jeA T 

on a single element T, where {V>j}jeA is the standard nodal basis of Vj-, we easily see that 

/ (AVwq-)-vds = / u>Tj(AVipj) • v ds 

JdfiinT £r^ T Jan^T 

= Y (wrj -wn) / (AVV>j) • v 
jeA T \{i} Jd ^nT 



ds . 



In the next step the matrix A is approximated by a piecewise constant matrix A7-. Summing 
up over all elements T which lie in the support of ipi, we have the relation 

/ {Aq-Vwq-) • Vipi dx = ( w Ti — w Tj) / (-A-tV^j) • vds . 
With the definition 

r d 

f-ij '■= 



■ ■ f 

-f- / (AV^)-uds, rriij >0, 

m ij JdUi 

, = 0, 



it follows that 



/ (AVu> r ) • Vipi dx « / (A r V-u; r ) • Vipi dx = V] ^ (w Ti - w Tj ) ^ . 
Jn Jn j£Az d^ 

Unfortunately, it is wellknown that in the case d = 3 the right equality does not hold. 
Nevertheless, the right-hand side - together with the above definition f|12|) of [iij - is senseful 
for d = 3, and thus this formula can be used for discretization. 

Consequently, in order to obtain a discretization for the case of a matrix- valued diffusion 
coefficient, it is sufficient to replace in the forms aj- and b-f the corresponding values of /ijj 
according to formula (]12p . 

REMARK 2 The really critical point in the discretization of diffusion- convection equations 
with matrix-valued diffusion coefficients consists in the appropriate choice of the stabilization 
mechanism in the situation where the eigenvalues of A are widely spreaded (cf. [AngOO], 
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2.3 The case of Donald diagrams 

Let us now consider a family of admissible (in the sense of FEM, cf. |Cia78t Ch. 2]) triangu- 
lations J- = {T}- Then, for any T E T with local vertices Zj = Xj. , ij £ At, j £ [l,d + 1]n, 
we define 

fi£ T := {x G T : (V& G [1, d + 1] N \ {j} : A fe (x) < A,(x))} , 



where Xj(x) is the j-th barycentric coordinate of x w.r.t. T. Define for i G A the sets 



fif := int |J nP T 

\T: dTBXi 

In this way, we get a family of Donald diagrams. 

Although it is possible to introduce a discretization like (jlip . we use the following version: 

a T (w T ; v T ) = (A(',w r )Vro r , Vu r ) + M^T! «r) + d T {w T - v T ), (13) 

where the forms 67-, dj- are defined analogously to ([9|). ([T0|) . In particular, 7^ G R is an 
approximation to (u ■ b)(-,wj-)\nj- 

In the case of a matrix- valued diffusion coefficient, we define nij analogously to (|12p 
but use it only in 67- to ensure a certain stabilization. The form aj- remains as it is, i.e. 
c$j-{wt;vt) ■= (A(-, wj-)Vwt, Vi>t)- 

3 Stability and a priori error estimates 

3.1 The case of a linear equation with a scalar diffusion coefficient 

In this section we give a short review of some wellknown properties of the schemes ([6]) and 
()13p for the case of a linear equation with a scalar diffusion coefficient. We start with the 
formulation of conditions with respect to the approximations mj and jij. 

(A2.1) is an approximation of the term J r Ads satisfying the following conditions: 

(i) < fj,ij < \\A\\i )0O) n, 

(ii) = fiji, 

(hi) /Ujj — m" 1 J r Ads < ChT\A\i jOC> n, where T is one of the simplices having the 
vertices X{,Xj, and C > is a constant independent of a,hx,i,j- 

(A2.2) jij is an approximation of the term m" 1 J T v-hds satisfying the following conditions: 

(i) < ||b||i )0 o,n, 

(ii) jij = 

(hi) jij — m^ 1 J r X u ij ' b)ds < C7iT|b|i |0O ,fi> where T is one of the simplices having 
the vertices Xj, Xj, and C > is a constant independent of b, hr, 

The subsequent results are extensions of the theory developed in |Ang91| , |Ang95b| , see also 
[KAH31 Ch. 6]. 
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THEOREM 1 (Discrete coercivity) Let a family T = {T} of triangulations be given, 
where in the special case of Voronoi diagrams (i.e. E = V) all elements T are self-centered 
and in the special case of Donald diagrams (i.e. H = D) the family is shape-regular. Moreover, 
let the assumptions (A2.1), (A2.2) be satisfied. Then, for ho > sufficiently small there exist 
two constants ao > and Hi > independent of h such that for all h G (0, ho] and vj- G Vj- 
the relation 

a T (v T ;vf) > aol^rll + ailbrllr 

holds. 

The a priori error estimate is based on this stability property and on the following con- 
sistency result. 

LEMMA 1 (Discrete consistency) Let a shape-regular family J- of triangulations {T} be 
given, where in the special case of Voronoi diagrams (i.e. E = V) all elements T are self- 
centered, and let the assumptions (A2.1), (A2.2) be satisfied. Then, if ho > is sufficiently 
small, for any element w G H V and any element vj- G Vr the estimate 

\a^-(L-j-w; v-r) — (— V • (AVw) + b • Vw + cw, LfVf)\ 

< C7i|H| 2 ,2,n [\vt\s + \\vt\\t] 

holds for all h G (0, ho], where C > is a constant which does not depend on h. 

The proof of the following theorem is a modification of the standard proof of Strang's first 
lemma. 

THEOREM 2 (A priori error estimate) Let a shape-regular family T of triangulations 
{T} be given, where in the special case of Voronoi diagrams (i.e. S = V) all elements T are 
self- centered, let the assumptions (A2.1), (A2.2) be satisfied and suppose that the solution 
u G V of problem (CP additionally belongs to W%(p). 
Then, for sufficiently small ho > the estimate 

\\u - u r \\s < Ch [|M| 2 ,2,fi + \ f\i, q ,n] 
holds for all h G (0, ho], where the constant C > is independent of h. 

3.2 The quasilinear case 

Due to the possible structural diversity of the nonlinearities in (pQ), in the nonlinear situation 
there is not such a relatively canonical theory as in the linear case. 

We mention here only a few papers which are concerned with the investigation of node- 
centered finite volume methods for nonlinear elliptic (or parabolic) equations and refer to the 
literature cited therein: [FLOl] . [CL05] . |EFG06j . 

4 A posteriori error estimates for nonlinear problems 

In this section we present the general approach that does not depend on the particular dis- 
cretization. 

The nonlinear primal problem we are interested in is given by 

u<EV: a(u;v) + a 5 (u;v) = (f,v) Vv G V. (14) 
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It represents the weak formulation of the originally given (accurate) boundary-value problem 
for a partial differential equation in a real Hilbert space V, where / is a linear functional on 
V and (/, v) denotes the value of / at the element v £ V. The forms a : yxl^->M and 
ag : V x V — >• R are linear in the second argument but may be nonlinear in the first one. In 
the context of the boundary- value problem (jl]), the left-hand side of (J3|) is written in (I14p as 
the sum a + a$, where a stands for a certain simplified problem and as represents a part of the 
equation which is to be neglected in the practical computations. That is, the discretization 
applies only to a in ()14p . The goal is to estimate the influence of both neglecting as and 
discretizing a and / with respect to a given output functional j : V — > R. 

EXAMPLE 7 Consider (CP with 

A(x, w) := e(x)\w\' y ( x >I, b(x, w) := bo(a:)M > c(x,w):=c (x), 

where e, cq, /, 7 : Q, — > R and bo : SI — > R d are smooth functions (satisfying certain additional 
conditions, in particular — 1 < 7_ < 7(2;) < 7+ < 00 on Q, for some constants 7-, 7+ )■ Then, 
for some constant elements wo,7o G R, we can set 

a(w;v) := (e\wo\ 10 Vw,Vv) + (|u>o| 70 t>o ■ Vw,J)) + (cow,v), 

and as(w;v) is the canonical error term with respect to the correct weak formulation of |7p. 

The directional derivatives of a(u; •) and as{u; ■) in u will be denoted by a'(u; •, •) and a' s (u; •, •), 
respectively. The form 

a'(u; w, v) := lim - [a(u + ew; v) — a(u; v)] 

is linear in w and v. The second and third directional derivatives are denoted by a"(u; •, •, •) 
and a"' (u; •), respectively. In the general case of a nonlinear output functional j, the 
corresponding dual problem we will use in the analysis is the following: 

z G V : a'(u;w,z) + a's(u;w,z)=j'(u;w) Vw £ V. (15) 

The solution z € V of the dual problem is called influence function for the particular choice 
of j ([AOOOj]). The primal solution u rn G V and the dual solution z m £ V of the reduced 
problems are given by 

u m £ V : a(u m ;u) = (/,t>) \/v £ V, (16) 
z m £V : a(u m ;w,z m ) = j'{u m ;w) Mw£V. (17) 

These variational problems will be formulated in terms of optimization problems. The primal 
and dual solutions will be expressed by the variables x := (u, z) £ X := V x V and x m := 
(u m ,z m ) £ X. In the variational space X, we consider the functionals 

L(x) := L m {x) + L s {x), (18) 
L m {x) := j(u) + (f,z) - a(u;z), (19) 
L s {x) := -a s {u-z). (20) 

The derivative of L applied to a test function y = (w, v) £ X is 

L'{x; y) = j'(u; w) - a'(u; w, z) - a' s (u; w, z) + (/, v) - a(u; v) - a s (u; v). 
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Obviously, the original primal and dual problems (|14p and (|15p and the reduced primal 



and dual problems (fl"6|) and (fTTj) consist of finding the stationary points x = (u, z) and 
x m = (u m , z m ) of L and L m , respectively: 

xeX: L'(x;y) = Vj/Sl, (21) 

x m eX: L' m (x m ;y) = Vj/GX (22) 



Furthermore, the target quantities are given by evaluation of L and L m at the following 
stationary points: 

j(u) = L(x), j(um) = L m (x m ). 

In order to balance the model and discretization errors, we have to include the discretization 
error in the analysis. To do this, let Vj- C V be a finite-dimensional subspace. Typically Vj- 
is a finite element space with respect to a partition T of the computational domain Q C R , 
d G {2,3}, where possible homogeneous Dirichlet boundary conditions are already included 
in the choice of the spaces V and Vf- Let aj- : Vj- x Vj- — > R be a nonlinear form which is 
different, in general, from the simple restriction of a to V7- x V7-, and denote by fj- : Vqr — > M. 
a linear functional which not necessarily coincides with f\v r - For instance, aj- and fj- may 
result from the finite volume discretization of a, / in (|14p according to Section [2J 
Then uq- m 6 V7- is the discrete solution of the problem 

u rm g V7- : a T (uTm;v) = (fr,v) VveV T (23) 

involving both types of error. The operators L and L m are still given by (I18p - (]20p . The 
difference lies in the definition of the discrete solution xj m = (u-j- m , zq- m ) £ Xj- := Vj x Vjr, 
where now Uf m satisfies (123H and Zf m is the solution of the following dual problem: 

z Tm G Vf : a'(u Tm ; w, z Tm ) = j'(u Tm ; w) Vw £ Vj. (24) 

In such a setting, the relations a(u-j- m ; v) = (/, v) and L' m {x-j m \ y) = are no longer valid 
for all v £ Vj resp. y G X7-. 

The target quantities are given by the evaluation of L and Lf m , where 

L Tm (x) := j(u) + (f T , z) - a T (u; z), (25) 

at the following stationary points: 

j(u) = L(x), j{u m ) = L Tm (x m ). (26) 

For the formulation of the error representation, we use the following notation for the primal 
and dual residual with respect to the reduced model and for test functions (w,v) G X: 

g(u Tm ;v) := (f,v)-a(u rm ;v), 

Q*(uTm; ZTm, w) := j'iUTm, w) - a'(UTm, W, ZTm)- 
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THEOREM 3 If a(u\-), ag(u;-) and the functional j(u) are sufficiently differentiable with 
respect to u, then we have 

j( u )-j( u Tm) = -asiuTm'jZTm) 

j ZTm ) - {fr, z Tm) ~ a(u Tm ; zrm) + ar{ u Tm\ ZTm) 

1 

+ 2 \Q\ u Tm\ z ~ irz) + Q*{uTm] ZTm, u - i-j-u)] 

&u > ZTm ) 

1 1 

- g a u Tm; ^r™ - *t*) ~ 2 ^' 

where e := (e u , e 2 ) := (u — tt7- m , z — z-r m ), W '■ V ~ > *s an interpolation operator, and the 
remainder R is given by 



R 



:= / a(l — a)L"' (x-j- m + o~e; e,e,e) da. 
Jo 



Proof: By ([26} , 



- j(u T m) = L(x) - L T m(xTm) 

= L{x) - L m {x Tm ) + L m {x T m) - L Tm {xTm) 

= L(x) - Lm(xTm) 

j ZTm ) - (/t, z Tm) - a(uTm] ZTm) + ar(uTm', z Tm), 

where the last step is a consequence of the definitions (f20l) . (|25j) . 

The first difference can be estimated as in the proof of [BE03, Thm. 2.1]: 

L(x) - L m (x Tm ) = L(x) - L(x Tm ) + L s (x Tm ) 

= / L'(x T m + 0-(x - X Tm );X - X Tm )do- + L S (x T m) 

Jo 

= - [L'(x Tm ; e) + L'(x; e) - R] - a 5 (u T m; z Tm ) 
with the above given remainder R of the trapezoidal rule. Since L'(x; e) = by (|2Tj) . we get 
L(x) - Lm(xTm) = -as(u Tm ; z Tm ) + ^ [L'(x Tm ; e) - R] . 

Furthermore, 

L'(x Tm ; e) = j'{uTm\ e u ) - a' {uTm\ e u , z rm ) - a' s (uT m ; e u , z Tm ) 
+ (f, e z ) - a(u Tm ; e z ) - a 5 (u Tm ; e z ) 
= 6*(uTm,;zTm,e u ) - a' s (uT m ;e u ,ZTm) + e(uTm',e z ) - as(uTm]e z ). 

Since the Galerkin orthogonality is violated, in general, we cannot use the standard argument 

= g(u Tm ; ZTm) = 0(uTm] itz) 
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to replace Zf m by i-yz in the third term. Here we can only make use of an analogous property 
of the dual problem ([23]) . i.e. 

= Q*(u Tm ;z rm ,UTm) = 6*(uTm,;ZTm,iru)- 

(Of course, if the dual problem is approximated by a finite volume method, too, then we have 
to argue as for the primal problem.) Thus we arrive at 

L'(x Tm ;e) = Q*(ur m ;z Tm ,u-iru) + g(urm;z-i T z) (27) 
- q( u t mj ZTm ~ ir z ) ~ a&{u Tm ; e z ) - a s (ur m ', &u, z Tm ). 

This gives the assertion. < 

In order to use numerically the error representation derived in Theorem (3] we have to 
approximate various terms. In particular, we will neglect the higher-order terms in e, namely 
the remainder R and the terms as(u-j- m ; e z ), a' 5 {uf m ; e u , Zf m ), cf. the related discussion in 
[BE03]. Furthermore, we have to approximate the interpolation errors u—ij-u and z — i-yz. An 
efficient possibility for doing this is the recovery process of the computed quantities by patch- 
wise higher-order interpolation expressed via the operator it- : Vqr — > formally, where 
Vt is a richer discrete space than Vj (see [BR01[ Sect. 5], [Ran05] Sect. 3.2]). For instance, 
in the case of triangles (d = 2) or tetrahedra (d = 3) and when Vj- consists of piecewise 
linear elements, quadratic interpolation may be used. For quadrilaterals and piecewise d- 
linear elements, the interpolation can be done on (i-quadratic elements. In order to preserve a 
sufficient high accuracy of the interpolation procedure, special care on elements with hanging 
nodes is required. 

The interpolation errors will be numerically approximated by 

z-i r z « i\zTm ~ ZTm, (28) 

U - %rU W l^j-UTm - UTm- 

Without the modeling error and in the case of conforming methods, this approximation is 
usually observed to be accurate enough. 

Taking into account that the residual (f(uf m \ z-j- m , v) vanishes with respect to a discrete 
test function v £ Vj~, we obtain from Theorem [3] the following approximate estimator con- 
sisting of three indicators: 

j(u) - j(u Tm ) » VT + Vm + Vnc, 

VT '■= 2 [s( U Tm,; ijZTm - ZTm) + Q*(ur m ; ZTm, l^-UTm)] , (29) 
; ZTm ) > 

Vnc ■= (f, ZTm) ~ {fT, z Tm) ~ a(uTm] ZTm) + ar{ u Tm\ ZTm)- 

The indicator rjr of the approximate estimator can be considered as the conforming contri- 
bution of the discretization, and the indicator r\ m measures the influence of the model. For 
complex models, the evaluation of rj m may be expensive. Often in practice the decomposition 
a + as is changed successively in such a way that portions of as are (locally) shifted to a. 
The indicator r] nc results from the nonconformity of the discretization method caused by the 
violation of the Galerkin orthogonality. The practical treatment of r] nc will be discussed in 
Section [5j 
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REMARK 3 According to t28\) , there are two ways for the treatment of the term 

g(u Tm ; z - i T z) - g(u Tm ; z Tm - Hz) 

occuring in Theorem^ Either we write it as g(wr m ;z — zf m ) (i.e. we reverse the splitting 
used in pT7|)j and replace then z by i^-z-r m , or we replace z — iq-z by i^Zf m — z-r m and ijz 
by z-r m . In both cases, we arrive at the same result: 

g{ur m ; z - ijz) - g{ur m ; zTm - irz) ~ g(ur m ; ij-zrm - zj m ). 

In order to use the information (I29p for changing locally the model or the discretization 
parameters (e.g. the mesh size), we have to localize the indicators. After that, an adaptive 
process has to be designed in order to balance the error sources. 

Regarding the localization of rjj- and r/ m , so here there are no new aspects. We refer, for 
instance, to [BE03]. 



5 Application to the finite volume method 

In the papers Ang91 , Ang92 , an extension of Babuska&Rheinboldt's a posteriori error 



estimates for finite element methods ( }BR78j ) to finite volume methods for linear diffusion- 
convection eqations has been proposed. In a subsequent paper (|Ang95a|), for a singularly 
perturbed model problem a modification was introduced with the aim to get two-sided bounds 
of the error such that the constants occuring in these bounds are independent of the pertur- 
bation parameter. In [AKT98] and [Thi99], residual type error estimates for finite volume 
discretizations of more complicated problems in two and three space dimensions have been 
presented. A rather general framework for the a posteriori estimation in various finite volume 
methods can be found in [V0I1O8J, however this paper is restricted to linear problems and es- 
timates w.r.t. the energy norm. In AnglO , dual-weighted residual error estimators for finite 



volume discretizations of linear diffusion-convection eqations have been described. Here we 
apply the results of the previous section to the nonlinear diffusion-convection problem. As 
a result, we get a posteriori estimates for errors of functionals depending nonlinear ly on the 
solution and for possible modeling errors. 

Interpreting a-r and fj- as the finite volume discretizations ([6]) of the forms a and / in 
(fl~4"|) , we first observe that the the estimators 777- and r? m depend only on the computed discrete 
solution but not directly on the structure of aj- and fx- Therefore, these estimators can be 
treated as in the (conforming) finite element case and we concentrate on the estimator fj nc . 
To simplify the presentation, we will write xr = (^Tj z t) instead of xrm = {uj 

Then, by definition, we have that 

(f,z T )-(fr,z T ) = ^{(/,z t )t-(/,^,t} (30) 

TeT 



{(f,ZT)T~ ^2 f^Ti 171 !}^ 



E 

TeT ieA T 



where (/, zj)t '■= / fzfdx and mj := meas^ (Qi n T) . Analogously, with 

(it,t{ u t; zt) 

■= ^2 z Til E ( 1 ~ ( UTi ~ u Tj)\ mfj + CiU Ti mJ 
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and cit(u-j-; z-r) resulting from the restriction of all integrals occuring in the expression for 
a(u-j-; zj-) to the domain of integration T, we have that 

a T (u T ; z T ) - a{u r ; z T ) = ^ {a TiT (u r ; z T ) - a T (u T ; z T )} . (31) 

TeT 

Putting (|30|) and ([3T]) together, we conclude that 

r]nc = 52 {(/) z t)t ~ 52 fi z Ti m f] + 52 i a T,T(ur; z t) ~ a T (u T ; z r )} . 
TeT i£A T TeT 

This is the starting point for the practical, localized computation. 

6 Analysis of the nonconformity estimator 

In this section we show for the case of a scalar diffusion coefficient A = AI with A: £1 x R — > 
R that r\ nc is order- consistent with the a priori error estimate (Theorem [2]). A precise 
formulation of this property is given at the end of the section. 
Using (p0|) . ([3T]) . we get the following decomposition: 



52 / if Z T ~ fiZTiW 



ieA 

+ 52 ZTi 52 My ("Ti - ^i)-^ 1 - {AVur, Vzr) 

ieA jeA; tj 
+ 52 52 ^ ~ r ijhij(uTj - UTijZTimij ~ (b • Vu r , Zr) 

ieA j'GA,; 

+ 52 / [<Hv>TiZTi - cur&rldx 

JVli 

ieA 

52 1 / ( zr ~ z Ti)^ + / (/ - fi)zn.dx \ 

■ -r U^i J 

+ 52 ZTi 52 ^ij( w Ti - u Tj)-j^- - (AVu T , Vz T ) 

ieA ieA, 1/1,1 
-52/ ( b ' ^ u t)( z T - ZTi)dx 

*t{i - / — } 



ieA 

#0 + #1 + ^2 + <>3 
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with 



5o := ZTi E ^( u Ti ~ UTj)-^ 1 ~ (AVu T , Vz T ), 

ieA jeA, tj 



51 := V] / if -b-Vu T -cu T ](z T - z Ti )dx, 

ieA 

5 2 := ^2 z Ti{ if ~ fi + (V • b - c)u T + CiU T i]dx - ^ "7 }• 

ieA jeA, 

^ 3 := EE/ [( r y' u Ti + (1 - nj)uTj)lij ~ {fij ■ h)u T ]z Ti ds. 



ieA jeA; 1 « 



Here we have used that b • V-ur = V • (hu-j-) — (V • b)u7". 
In the case of Donald diagrams, <5o = 0. 

In order to treat do in the case of Voronoi diagrams, we introduce a piecewise constant 

(w.r.t. T) approximation Aj to A by An-W '■= : — r / Adx, T G T. Then we can write 

meas d (T) 7 T 



E z ^ E f ^' ~~ [ Ards ) ( UTi ~ ^i)^r 
ieA jeA, V Jt * J " : ' 

+ E ^ E ( I A rds) UTl ~ UT3 - (AVu T , Vz r ). 



ieA jeA 

It is wellknown that, for arbitrary uj-, zj- G V7-, 



E ^ E ( / ^1 = (4rv«r, v ^)- 

*eA .;. \. V 7 ' / 



Hence 



*° = E ZTi E ( ^ ~~ Ards ) ( uri ~ ut ^~i 1 + ^ Ar ~ A ) Vur > Vzr )- 



ieA jeAi 



Since both Vu-f, Vzj- are constant on every element T G T, the second term vanishes. By a 
symmetry argument, we arrive at 

60 = \ E E f ^' / Ards J ~ u Tj)( z Ti - ZTj) 1 ^ 1 - 

ieA jeA, V "'" - /r / 

Now the Cauchy-Schwarz inequality implies 

N < UeE^-^-/ ^ (^-n Tj ) 2 ^} 

2 [ ieAieAi V m ^ ^ / ^ J 



x < 



1/2 



EE(^-^)« 

ieA jeA, lJ 



L. Angermann: Error estimates for functionals. March 3, 2013 



17 




fig n fi, 



fig' n fii 



Figure 2: The auxiliary simplices in the case d = 2 for the Voronoi diagram 



The last factor can be bounded by C\\zq-\i t 2,n-, therefore we get 



< Cmo\zr\i,2,n, 



(32) 



where 



2 J2 



with 77^ := 7 ( /% - — / A T ds J 



2 Tly 



Setting g := f — b ■ Vtif — raj- and <5ii := 5(27- — z^i)dx, we can write (cf. Figure [2] for 
the case d = 2): 



hi =^2 ^2 T 9 ( ZT ~ z Ti)dx. 



je^iT&T:mf,>0 V 



On each simplex T, it holds 



Z T = Z T j + Vz T ■ (X - Xi), 



where Vzj- is constant on ttj-. 
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It follows that 



hi 



E E T 9 ^ zr ■ ( x ~ Xi "> dx 



< 



E E 

jeA* TeT:mT.>0 Jn Ii nQ 



\g\\\V z T \\\\x - Xi\\dx 



< < 



1/2 



E E 



| | 2 1 1 a? — Xj|| 2 dx 



<< E f T \^ z r\\ 2 dx\ 

\ lJ ) 

)l/2 



1/2 



prli,2,Qi 



< vii\zr\i,2,ni, 



where 



Thus we arrive at 



Vi'-^Vii with rfii := ^ ^ /i| / £> 2 cte. 
#i < r/i|^7-|i,2,n- 

For the third term 82, with 

0» := / [/ - fi + (V ■ b - c)u r + Qn Ti ]dx - V] "7 •"<./'"/./• 
Jtti 771. 



we have 



Because of 



(33) 



$2 = ^2 ZtA- 

ieA 

ZT'fii < »72i I ^Ti I \/"*i> 
„2 ^2 



where r\2i := it follows with 772 := X^eA 7 ?^ that 

^2 < ^Ikrllr- 

In view of the equivalence of the L2-norm and the lumped L2-norm on Vqr, we obtain 

82 < C 2 r]2\\zr\\o,2,n- (34) 
For the remaining term 83 we have, by a symmetry argument, that 

iGA 
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where 



with 



S M '■= 2 E / ^Tij{zn ~ z T j)d& 



Crij ■= [rijU Ti + (1 - rij)uTj]%3 - {vij • b)u T . 



In view of zn — zq-j = dijiy^ ■ Vzf) on we get 



j/'gAi TGT :m^>0 Fi J 

"J 

It follows (remember that z/^ • Vz-j- is constant on rg and Vzj- is constant on 

fig n Hi) 



<5 3 * < 



2 E *i E 

ieAi TeT:mJ.>0 

?E E 



(Tijds 



\Vz T \ 



jeAi TeT ■. >o A / measrf ( fig n fi, 



»3 



By Cauchy's inequality, we have 



hi < - < 



1/2 



4 



EE-, . . 

i6Ai TeT ■■mT->Q meaSrf ( fi « R °* ) V Jr 5 



(Tijds 



PT i,2,a 



< %ikr|i,2,ni 



where 



^1 



4 



JE E 

4 ^ ^ m£ \ /r' 



jeA, TeT:m T >0 meas d ( ^« n ) \ Jr 5 



(Tijds 



(Tijds 



Thus it holds that 

#3 < %|^r|i,2,c- 
Summarizing the estimates ([32]) - ([35]) . we obtain 

Vnc < (CiTjo + r]i + %)kr|i,2,n + C , 2 r?2lkrllo,2,n, 

where the indicators have the following structure: 

( \ i/a 



»7! 



E 

. ieA 



i€ {0,1,2,3}, 
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where 



VOi 



Vu 



V2i 



V3i 



1 

- < 

2 



X) ( ^ - — / A rds ) (u Ti - UTj) 2 ^ 11 > 



1/2 



in case of Voronoi diagrams and ??oi = in case of Donald diagrams, 

( \ 1/2 



E h T I T [f-h-Vu T -cu T ] 2 dx\ , 

j'GAi TgT: rnjj>0 ■ /n 5 n °< j 

/ [/ - /« + (V • b - c)u r + auri]dx - ^ UTHiji 



5 X X ~T\ f t ^ TijUTi + ( X ~ r v) u Tjhij ~ {vij ■ h)u T ]ds j 



1/2 



REMARK 4 (%) We mention that all the indicators rji can be rewritten in such a way that 
the resulting local indicators are related to the elements T E T ■ 

(ii) It can be shown that the indicators r\i are order-consistent with the a priori error 
estimate (Theorem^) in the following sense : 

If f £ Wg(Q.) with some q > d and u G W 2 2 (^)i then there is a constant C c > such that 

3 

$><c c MI|u|| 2 ,2 + ||/M, 

1=0 

see \Ang92 , Thm. 4] for a special case. 



7 Conclusions and perspectives 

We derived an estimator for measuring simultaneously two types of errors, modeling and 
discretization errors, with respect to user-defined output functionals. The approach is formu- 
lated for stationary nonlinear partial differential equations involving complex models. The 
main focus was on the consideration of discretization methods which do not possess the prop- 
erty of Galerkin orthogonality. For the example of node-centered finite volume methods, by 
localization of the estimators we presented local error indicators which allow for local mesh 
refinement and local model modification. 

In future work, the described framework will be the basis for a more detailed analysis of 
problems with tensor-valued diffusion coefficients and dominating convection and for algo- 
rithms which balance the indicators corresponding to the different sources of error. 
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